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Consider  a  quadratic  function  g  of  the  form 
g(x)  *  x^Cx  +  2d^x 

with  d  an  n- vector  and  C  a  synmetric  n  by  n  matrix. 

A  classical  manipulation  reduces  g  to  the  form 

g(x)  =  xT  C  x 

T 

where  x  =  (l,x)  and  C  =  (°  £  )  ,  syranetric.  We  may 
therefore,  without  loss  of  generality,  restrict  our  analy¬ 
sis  to  quadratic  functions  of  the  form: 

g(x)  =  xT  Cx  ; 


I 


in  fact,  C  need  not  even  be  assumed  symmetric  in  the 
following  developments. 

Note  that  no  assumptions  concerning  definiteness  of  g  is 
required  either. 

We  now  define  a  Polaroid  bifunction  f  which  polarizes  g  : 

f(x;y)  -  xT  Cy 

Property  1:  f  is  simple. 

Proof:  Imnediate  from  the  definition  1  since  f  is  bilinear. 

Q.E.D. 

Polaroid  sets  are  defined  (see  [1))  in  terms  of  a  polaroid  bifunction 
f  ,  for  an  arbitrary  (closed)  set  P  and  a  parameter  kcf(P;Bn)  .  For 
our  present  purposes,  we  only  need  a  definition  related  to  the  above 
bilinear  bifunction  f  : 

jc 

Definition  3:  The  simple  Polaroid  set  P  (k)  is  defined  by 
P  (k)  =  {y  1  f(x;y)  =  xTCy  <  k  ,  VxeP} 


3 


2)  Convexity  and  validity  of  the  simple  polar olds. 


Proof: 


* 

The  simple  polaroid  P  (k)  is  a  convex  set. 

O  1  1  o  ^ 

Let  y  =  Xy  +  (l->  y  ,  with  y  and  y  e  P  (k)  , 


3 

and  X  e  [0,1]  ;  we  need  show  y  e  P  (k)  ,  i.e.  VxeP 


xTCy3  =  XxTCyX  +  (i-X)  x^Cy2  <  Xk  +  (1-X)  k  =  k 

Q.E.D. 

* 

Furthermore,  the  polaroid  P  (k)  can  be  considered  as  a  valid  outer- 
domain.  in  the  sense  chat  it  contains  no  point  xcP  with  a  value  g(x) 
greater  than  k  ; 


Proof: 


3:  P  n  P  (k)  c  levk  g  . 

*  T 

Take  y  e  P  ("I  P  (k)  :  then  x  Cy  <  k  ,  VxeP  ;  in  particu¬ 
lar  for  y  e  P  one  has 


y  Cy  =  g(y)  <  k 


Q.E.D. 


Note  that  no  particular  assumptions  concerning  P  need  be  made; 
thus  the  cutting  planes  described  below  are  valid  for  any  subset  of  P 
as  well;  this  is  clearly  a  useful  property  in  the  context  of  integer 
quadratic  programming,  for  instance. 

2.1  Cutting  plane 3. 

Suppose  that  P  is  polyhedral  and  characterized  by  the  following 


x  -  b  -  At  >  0  ,  t  >  0 


explicit  format: 


the  n-'vector  b  and  the  n  by  n  matrix  A  are  derived  from  the  full 
linear  programming  tableau  corresponding  to  an  arbitrary  feasible  basic 
representation ,  by  choosing  the  rows  and  columns  relevant  to  the  initial 
variables  x  ;  t  denotes  the  n-vector  of  n on-basic  variables. 

Consider  now  the  mathematical  programming  problem 
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Since  it  is  obviously  desirable  to  generate  as  deep  a  cut  as  possible, 

ic  • 

one  will  determine  the  value  X^  >  0  ,  j«N  which  corresponds  to  the 
intersection  of  the  ray  u^  with  the  boundary  of  P  (k)  (i.e. 

f (xju*)  =  k)  ;  thus  one  requires 


•P  T  —  •fc  m 

x  Cu,  =  x  C(b  -  X,  a.)  =  k  ; 
j  J  J 

Assuming  thut  there  e  ists  a  point  u^(X)  ,  Xe  [0,1]  which  lies  in 
Int  P  (k)  ,  we  can  infer  that  the  intersection  point 

^  -j^  ^  ^ 

u^  =  uj(^j)  *  bd  P  (k)  with  Xj  >  0  is  imique.  (This  assumption  is 

-T 

easily  satisfied  by  choosing  k  >  b  C  b  ,  for  instance.) 

The  intersection  problem  merely  reduces  to  finding  one  point  xep  such  that 
<£  _  ^ 

x  C(b  -  Xj  a^)  =  k  ;  this  can  be  solved  by  the  following  parametric  linear 
programming  scheme: 


0: 


i 

Set  X.  =  X°=  min  - 

J  J  su  > 0  sa 


(simplex  criterion) 


Step  _1: 


(Any  other  x”  >  0  can  also  be  used  as  initial  value.) 
Solve  the  L.P.: 


T  - 

max  z(x)  =  x  C(b  -  X^  a^)  ,  X^  fixed 


subject  to  xcP 


Let  x,  z  =  z(x)  denote  the  optimal  solution. 


,  * 

and  set  Xj  =  Xj  +  AXj 

"ic  /v* 

If  Xj  =  Xj  (i.e.  AXj  =  0,  because  z=*lc)  then  stop! 

Otherwise  update  the  value  Xj  by  setting 

,new  *  ,  .  , 

Xj  :  *  Xj  ,  and  go  to  Step  1. 

It  can  be  shown  that  the  polyhedral  nature  of  the  set  P  implies 
finite  ness  of  the  above  iterative  procedure. 


! 
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2.2  Applications :  In  conclusion,  we  note  that  a  deep  polaroid  cut  can  be 

generated  for  the  optimization  of  a  quadratic  function  sub¬ 
ject  to  linear  constraints.  This  cut  can  be  implemented 
in  a  variety  of  circumstances,  for  various  types  of  problems  and  in  vari¬ 
ous  kinds  of  algorithms. 

-  Problem* ;  The  present  polaroid  cuts  apply  to  a  more  general  class  of 

problems  than  die  strictly  concave  quadratic  programs  of  [4] 
since  definiteness  is  not  required;  furthermore,  in  the  concave 
case,  they  can  be  shewn  (see  [2])  to  dominate  uniformly  Hoang  Tuy's  cuts. 
Moreover,  our  polaroid  cuts  apply  to  any  kind  of  problem  where  F 
represents  a  superset  of  the  feasible  set;  typical  applications  of  this 
type  occur  when  the  original  problem  contains  some  extraneous  conditions 
(Integrality,  zero-one,  extreme  point  conditions,  etc.  ...  see  [5,6]). 

-  Algorithms :  A  first  straightforward  implementation  can  be  thought  of 

in  terms  of  a  method  of  cuts  accompanied  by  a  dual  opera- 
tionsto  restore  primal  feasibility  (this  procedure  is  well- 
known,  particularly  in  integer  programming).  As  usual,  thx3  approach 
may  present  here  some  convergence  difficulties. 

Because  they  do  not  require  definitness,  the  present  cuts 
were  also  found  useful  in  curtailing  a  face  decomposition  tree  search 
(see  [7,8]).  Clearly  they  can  also  be  used  in  any  other  branch  and  bound 
approach. 

The  next  section  establishes  further  properties  of  the  simple  polaroids 
P  ,  and  relates  them  to  further  types  of  algorithms  (see  also  [3,9]). 


3)  The  facial  structure  of  the  simple  polaroids  P  . 


Proposition  2:  Suppose  C  non-singular;  then  to  every  point 

x  e  bd  (conv  P)  there  corresponds  a  hyperplane  (halfspace) 
xTCy  (g)  k 

which  supports  P  (k)  at  a  point  y  e  bd  P  (k)  . 


Proof: 


Proof: 


|f  ^ 

By  definition,  P  (k)  lies  in  the  half3pace  x  Cy  <  k  , 

for  any  given  xeP  .  Consider  now  a  hyperplane 
T 

x  u  (<p  k  supporting  the  set  conv  P  at  the  point  x  --  x  ; 
let  y*  =  C**^u  . 

T  ★  T  ★  * 

Then  x  Cy  =  k  and  x  Cy  <  k  VxeP  ;  thus  y  lies 

★  x  * 

on  the  boundary  of  P  (k)  and  x  Cy  <  k  supports  P  (k) 


at  y  =  y  . 


Q.E.O. 


Corollary  2.1:  Suppose  C  singular;  let  y  define  an  objective  function 

T  * 

z(x)  =  x  Cy  which  satisfies: 

max  z  —  z(x)  =  k  ,  (x*P)  . 

x«P 

~T  *  * 

Then  the  hyperplane  x  Cy  (<)  k  supports  P  (k)  at  y  =  y 


Similar  to  that  of  Proposition  2. 


Suppose  P  closed,  convex,  polyhedral;  i.e.. 


Q.E.D. 


P  =  {x  |  Ax  <  b]  ,  with  suitable  A  and  b  .  Let  a  face  F 
of  P  be  represented  by  the  system  ax  =  3  where  a  and 
B  are  obtained  from  A  and  b  respectively,  .y  deleting 
the  appropriate  rows.  Then  to  F  there  corresponds  a  face 
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'fc'ic 

F  of  P  (k)  whioh  is  characterized  by 

F*  =  {y  |  Cy  =  orTp,  ,  *  k  ,  \l  >  0} 

T  T 

Proof:  By  hypothesis  one  has  for  ^  >  0  :  n  ax  =  p,  0  ,  Vx«F  and 

T  T  T 

by  transposition  isa|i  =  P|i  ;  thus  if  ^  is  chosen  in 

X  * 

such  a  way  that  3  |x  =  k  ,  one  obtains  for  P  (k)  the 

T  T 

supporting  hyperplanc(s)  x  u  (<p  k  ,  with  u  =  a  p.  ,  p,  >  0  . 
And  for  all  y  such  that  Cy  —  u  one  has 

xTCy  fcg)  k  ,  VxeF  indicating  that  y  e  F*  c  bd  P*(k) 

Q.E.D. 

T 

Note  that  the  relation  Cy  =  a  p,  ,  defines  the  respective  dimensions 

* 

of  corresponding  faces  F  and  F  ; 

Corollary  2.3;  Let  C  be  non-singular.  Then  to  a  p-dimensional  face  F  of  P 

*  ■* 

there  corresponds  a  (n-p-1)  dimensional  face  F  of  F  (k)  . 

*  —1  x  Y 

Proof :  One  has  F  ={y  |  y  =  C  <*  p  ,  p  >  0  3  H  =  k}  ;  by  assump- 

-1  T 

tion  the  matrix  C  a  has  rank  (n-p),  but  the  (n-p)  vec¬ 
tor  p.  has  only  (n-p-1)  degrees  of  freedom  on  account  of 
3  (X  k  •  Q.E«D. 

(Note:  as  in  corollary  2.1,  a  similar  statement  can  be  made  here  when 

* 

C  is  singular,  indicating  that  the  dimension  of  F  is  >  (n-p-1)  ). 

An  interesting  example  of  corollary  2.3  is  obtained  for  p  =  n  -  1  , 
i.e.  the  facets  of  P  ,  which  correspond  to  the  rows  of  the  matrix  A 

(defining  inequalities);  indeed,  one  obtains  the  (n-n+1-1)  =  0  dimensional 

*  *  * 
faces  F  ,  i.e.  the  vertices  of  P  (k)  .  Hence  the  vertices  of  P  (k)  can 

be  insnediately  given  by  the  relation 


,Tjr  -»  lCr-> 
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v  =  (c”*arT)  r~  where  cr  is  any  (single)  row  of  A  , 

q  bq  q 

This  property  holds  in  particular  for  the  classical  polar  sets  where 
the  matrix  C  is  assumed  definite  positive  (it  is  the  identity  matrix  for 
instance  in  [10]).  It  is  used  in  [3]  to  construct  a  new  algorithm  for 
strictly  concave  quadratic  programing,  which  is  based  on  the  vertices  of 
the  polar  set. 

In  this  context,  the  results  of  section  3  can  be  seen  to  introduce  a  new 


tool  for  the  following  problems: 


•  semi-,  and  indefinite  quadratic  programs ;  a  linear  term  may  produce 
a  singular  matrix  C  ;  in  an  algorithm  based  on  extreme  points  v  (like 
[3]),  one  may  replace  here  the  vertex  v  by  any  point  w  satisfying 


k_ 

b 

q 


This  holds  for  an  arbitrary  indefinite  quadratic  problem. 

*  OQP~lin**r  extreme  point  and/or  integer  programming:  when  the  feasible 
set  S  ia  only  a  subset  of  P  (of  integers, for  instance),  we  know 

“it 

that  P  (k)  is  still  valid  because  the  inclusion  theorem  (see  [1]): 

*  * 

S  c  p  =»  s  3  P 

may  be  invoked;  simple  polaroids  may  therefore  be  expected 
to  contribute  to  the  solution  of  such  less  tractable  mathematical 
programming  problems. 

In  conclusion,  we  establish  a  property  which  helps  visualizing  the  cor- 
* 

respondence  F  -♦  P  which  turns  out  to  be  involutory  under  the  proper 
assumptions : 

Property  4:  Assume  C  symmetric  and  non-singular;  then 

(P*(k))*(k)  =  conv  ({0}  U  P) 


(?  is  assumed  closed) 


Proof: 


First  let  us  note  that  0  *  P  (k)  ;  furthermore,  we  know 


that  a  polaroid  is  a  closed  convex  set  (property  2);  hence, 

"Me 

P  must  be  convex  and  contain  the  origin. 

Let  Q  =  conv  ({0)  U  P)  ;  we  show  that  Q*(k)  =  P*(k)  : 
the  inclusion  theorem  (see  [1])  implies  from  P  c  Q  that 

it  it  ^ 

Q  c  p  ;  furthermore,  for  y  e  P  (k)  one  has 
T 

x  Cy  <  k  ,  VxsP  .  An  arbitrary  point  q  e  Q  can  be  des¬ 


cribed  by 


q  -  Xx  +  (1-X)0  =  Xx  ,  X  e  [0,1] 


hence  q^Cy  =  Xx^Cy  <  xTCy  <  k  ,  VqeQ 

^  it 

this  implies  that,  for  any  y  e  P  (k)  „  one  has  y  e  Q  (k); 

*  * 

thus  P  —  Q  . 

it  it 

It  now  remains  to  show  that  (Q  (k))  (k)  ~  Q  , 

*  * 

By  definition  the  convex  set  (Q  (k))  (k)  is  generated 
by  the  following  family  of  hyperplanes  (half3paces) 
uTx  (<p  k 

where  uT  =  yTC  (i.e.  u  =  CTy)  .  v$*h  y  e  Q*(k)  =  P*(k)  . 

Consider  now  a  family  of  supporting  hyperplanes  (half- 
X 

spaces)  v  q  (<)  k  which  defines  Q  ;  the  existence  of 
such  a  family  preseuts  no  difficulty  since  every  convex 
set  can  be  considered  to  be  the  intersection  of  a  collection 
of  halfspaces.  On  the  other  hand,  one  has 
qTCy  <  k  ,  VqcQ  ,  y  e  Q*  (k)  ; 

since  q^Cy  =  y^C^q  =  v^q  <  k  ,  we  note  that  the  vector 


v  *  Cy  ,  y  c  Q  (k)  ,  generates  a  hyperplane  in  the 
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m 


lattes 


family;  furthermore,  the  definition  of 


Q  (k)  indicates  that  every  v  of  this  family  generates 
an  element  y  =  C  \  t  Q*(k)  which  is  uniquely  defined 

when  C  is  non-singular;  hence,  assuming  symmetry 

X  ** 

(i.e.  C  =  C  )  ,  one  sees  that  both  sets  Q  and  Q 


are  characterized  by  the  same  family  of  hyperplanes  (half- 

X  * 

spaces)  u  x  (^)  k  ,  with  u  =  Cy  =  v  ,  y  e  Q  (k) 

**  .  Q.G.D. 


implying  Q  =  Q 


iYra«. 
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